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Abstract 

By means of the Drinfeld twists, we derive the determinant representations of the 
partition functions for the g/(l|l) and gl(2\l) supersymmetric vertex models with do- 
main wall boundary conditions. In the homogenous limit, these determinants degen- 
erate to simple functions. 
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1 Introduction 



The six vertex model on a finite square lattice with the so-called domain wall (DW) boundary 
condition was first proposed by Korepin in pQ, where recursion relations of the partition 
functions of the were derived. In [2], [3] it was found that the partition functions of the DW 
model can be represented as determinants. Taking the homogenous limit of the spectral 
parameters, Sogo found that the partition function satisfies the Toda differential equations 
[I]. Then by using the equations Korepin and Zinn- Justin obtained the bulk free energy of the 
system [5]. The determinant representations of DW partition functions are uesful in solving 
some pure mathematical problems, such as the problem of alternating sign matrices |6j. By 
using the fusion method, Caradoc, Foda and Kitanine obtained the determinant expressions 
of the partition functions for the spin-fc/2 vertex models with the DW boundary conditions 
[7|. In [8], Bleher and Fokin obtained the large N asymptotics of the DW six- vertex model 
in the disordered phase. 

Supersymmetric integrable models based on superalgebras are important for describing 
strongly correlated electronic systems of high T c superconductivity [H QUI dH [321 EH EEH EES]- 
In this paper, we investigate the gl(l\l) and g/(2|l) supersymmetric vertex models on a 
N x N square lattice with DW boundary conditions. By using the approach of Drinfeld 
twists [T6l [T8l [19] , we derive the determinant representations of the DW partition functions 
for the two systems. We find that the partition functions degenerate to simple functions in 
the homogenous limit. 

The framework of the paper is as follows. In section 2, we obtain the determinant 
representation of the partition function of the g/(l|l) vertex model with DW boundary 
condition. Then we derive the homogenous limit of the partition function. In section 3, 
we solve the g/(2|l) supersymmetric vertex model with the DW boundary conditions, by 
deriving its DW partition functions exactly. In section 4, we present some discussions. 

2 3^(1 1 1) vertex model with DW boundary condition 

In this section, we study the DW boundary condition for the gl(l\l) vertex model onaAfxiV 
square lattice. 
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2.1 Description of the model 



Let V be the 2-dimensional irreducible g/(l|l)-module and R G End(V <E> V) the i?-matrix 
associated with this representation. In this paper, we choose the FB grading for V, i.e. 
[1] = 1, [2] = 0. The the i?-matrix satisfies the graded Yang-Baxter equation (GYBE) 
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(2.1) 



where Rij = Rij{\, Xj) with spectral parameters specified by Aj. Explicitly, 
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where 



a 12 = a(Ai, A 2 ) = 
C12 = c(Ai, A 2 ) 



Ai — A 2 
Ai - A 2 + 77' 
Ai - A 2 - 



612 = 6(Ai, A 2 ) = 



rj 



Ai - A 2 + 77' 



(2.3) 



Ai - A 2 + 77 

and 7/ is the crossing parameter which could be normalized to 1 in the rational cases we are 
considering. 

In the following, we study the (?/(l|l) 6- vertex model on a 2-d N x N square lattice. At 
any site of the lattice, there may be a vacuum or a single fermion p. A vertex configuration 
in the lattice is constructed by two nearest particle states in a horizontal line and two nearest 
particle states in a vertical line. For the present model, there are altogether 6 possible weights 
corresponding to a vertex configuration. The possible configurations and their corresponding 
Boltzmann weights Wi are given as follows: 
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Figure 1. Vertex configurations and their Boltzmann weights for the gl(l\l) vertex model. 
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Figure 2. DW boundary condition for the gZ(l|l) vertex model. 



One may check that these weights preserve the fermion numbers. If we assign horizontal 
lines parameters {Xj} and vertical lines parameters {zk}, and let 

Wi = l, w 2 = c(X,z), w 3 = w± = a(\,z), w 5 = wq = b(X, z), (2.4) 

then the Boltzmann weights correspond to the elements of the gl(l\l) i?-matrix (12. 2\\ . 

We now propose the g/(l|l) vertex model with DW boundary condition. Similar to [US], 
if the ends of the square lattice satisfy the special boundary condition, that is states of all 
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left and top ends are vacuum while states of all right and bottom ends are single fermion 
p, we then call the boundary condition the DW boundary condition. 

A configuration with the DW boundary condition is shown in figure 2. In this figure, Aj 
and Zj are parameters associated with the horizontal and vertical lines, respectively. The 
partition function for this system is then given by 

where the summation is over all possible configurations, rtj is the number of configurations 
with the Boltzmann weight By means of the R- matrix, the partition function may be 
rewritten as 

N-l N-l N N N N 

i=0 j=0 fc=l 1=1 j=l i=l 

Here pu\ (<f>u)) stands for the state p (0) in ith horizontal (vertical) line. 

Following Korepin, we use the transfer matrix to arrange the partition function. Define 
the monodromy matrix Tfc(Afc) along the horizontal lines 

^fc(Afe) = i?fc 5 Ar(Afc, ZN)Rk,N-i{^k, zn-i) ■ ■ ■ Rk,i{<^ki z i) = \ r<(\\ n(\ \ ■ (2-7) 



_C(X k ) D(X k ) 

therefore, in terms of T^(Afe) the partition function is equal to 

N-l N-l N N N 

i=0 j=0 k=l j=l i=l 

N-l N 



II (l| ( ^ i) C(A 1 )^(A 2 ) . . . C(X N ) II |0> (4) . (2i 



i=0 i=l 

Here we have introduced the notation 

|0>=(;), ID=(J) (2.9) 
to denote the particle states 0, p, respectively. 

2.2 Exact solution of the partition function 

In this subsection, we will compute the partition function (12.81) for the gr/(l|l) vertex model 
with DW boundary conditions by using the Drinfeld twist approach. 
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Let us remark that even though we believe that Drinfeld twists do exist for all algebras, 
so far only those related to A-type (super) algebras and to XYZ model have been explictly 
constructed pU QH [191 1201 EQ . 

Let a be any element of the permutation group Sjy. We then define the following lower- 
triangular matrix [191 I2Q] 

* N 

F h ... N = E E II p 53)) 5 k a ' a °) R °-N> (2- 10 ) 

(TG5jv Q! ct (1)— a CT (jv) j = l 

where P^ has the elements {Pn,\) m n = ^am^cm at the fcth space with root indices a = 1,2, 
the sum J^* is taken over all non-decreasing sequences of the labels a a (i)- 

Oia(i+i) > a<r(i), if c(« + 1) > cr(i), 

a<j(i+i) > Qf ff( i), if cr(i + 1) < er(i) (2.11) 
and the c-number function S^c, a, a a ) is given by 

S(c,a,a a )=expU ^ (l - (-l) [<Mfc)I ) ^ W ^ m ln ( X + M*Mo) } ■ ( 2 - 12 ) 

In (12.101) . R" 2 N is the iV-fold F-matrix which can be decomposed in terms of elementary 
F-matrices (I2.2p by the decomposition law 

-^L.JV = Ra'(l...N)Rl...N- (2.13) 

We showed in [18] that the F-matrix is nondegenerate and satisfies the relation 

F a (i...N)(z (T (i), . . . , z u ( N ))R[^ N (zi, . . . , z N ) = Fi...jv(zi, . . . , z N ). (2-14) 

The nondegeneracy of the F-matrix means that its column vectors form a complete basis, 
which is called the F-basis. The nondegeneracy also ensures that the F-matrix is invertible. 
The inverse is given by [19] 

/••,.iv vIl A '/ ( 2 - 15 ) 

i<j 

with 

Aij = diag ((1 + Qj)(l + Cji), dji, aij, 1) (2.16) 
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and 



jV 

&■))■ 



ct€<Sjv a CT (i)— a CT (jv) J—l 

(2.17) 

Here the sum Yl** is taken over all possible on which satisfies the following non-increasing 
constraints: 

< ®<r(i), if <?{i + 1) < ° r (0> 
a a (i + i) < a a (i), if cr(z + 1) > a(i). (2.18) 

Working in the F-basis, the entries of the monodromy matrix (12.7ft can be simplified to 
symmetric forms, e.g. the lower entry C(A) becomes [IH] 



C(A) = Fi2...nC{X)F 1 2 _ n 

N 

= E & (A,^)4)®^diag(2a(A,^),l) . ) , (2.19) 
i=i 

where i?^ (a, /3 = 1, 2; Z = 1, 2, . . . , N) are generators of the superalgebra g/(l|l) at the site 
/. From (12.191) . one can see that in the F-basis all compensating terms (polarization clouds) 
in the original expression of C(X) in terms of local generators disappear from C(A). 

Applying the F-matrix and its inverse to the states |0)(i) ® • • • ® |0}(jv) and (1|(at) <£> ■ ■ ■ <8> 
(l|(i), we have 

Fi... N \0)(i) ® • • • ® |0)(tv) = |0) ( i) ® • • • <g> |0) (7V) 

(1| ( jv) ® • • • ® (11(1)^7.^ = (1|(at) ® • • • ® (l|(i) n( 2a (^' z i)) _1 - ( 2 ' 2 °) 
Substituting (12.191) and (12.201) into the partition function (12. 8p . we obtain 

N-l N 

z N = n< i i^-o cr ( A i) c '( A 2)...c(Ajv)ni o >w 

i=0 i=l 
N-l N 

= IT (M(N-i)F^] N F 1 ^_ N C '(Ai) . . . C(Ajv)F 1 7 1 Ar Fi...jv |0)(j) 

i=0 i=l 
iV-1 TV 

= n( 2a ^> ^o)- 1 n <ii(iv-,)C'(A 1 ) . . . c(\ N ) n io> (0 

i<j i=0 i=l 



N-1 

JJ(2a(^,^))- 1 JX< 1 l(iv-i) 

i<j 



i=0 



N 



x2^- £ B N (Xi, . . . , Ajvl^ii, • • • > z i N )E(ii) ■ ■ ■ E \i n ) \\ |0)(i) 

i=i 



■tl<...<«]V 



j J a Zj)B N (\i, . . . , Ajvl^l, • • • , 2!jy), 



(2.21) 



i<j 



where 



N 



N 



=k+l 



B N (X 1 ,...,X N \z 1 ,...,z N ) = ^ sign(a)Y\b(\ k ,z a ( k )) J] a{X k ,z a{V) ] 

ctG<Sjv k=l 

= det£({A}, {z}) 

with £>({A}, {z}) being a N x N matrix with elements 



(2.22) 



a-l 



(B({X},{z})) af3 = b(X a ,zp)l[a(X„ze). 



(2.23) 



7=1 



Therefore the partition function of the g/(l|l) vertex model is given by the following deter- 
minant of the N x N matrix, 



Z N = Ha^iz^z^detBdXy^z}). 



(2.24) 



i<j 



2.3 Homogenous limit of the partition function 

In this subsection, we discuss the homogenous limit, i.e. when Ai = A 2 = . . . = A^ and 
Z\ — z 2 — . . . — Zjv, of the partition function (I2.24p . 

For later convenience, we rewrite the inhomogeneous partition function more explicitly 

as 



Z N = U 



i<j 



z, - Zj + rj 

Z{ Zj 



X 



b(X 1 ,zi) 
b(X 2 , zi)a(Xi_, zi) 

N-1 

&(Ajv, z\) Y[ a(A 7 , z{] 

7=1 



b(\i,z k ) 
b(X 2 , z k )a(Xt, z k ) 

N-1 

b(X N ,z k ) Y[ a(A 7 ,Zfc) 

7=1 



b(Xi,z N ) 
b(X 2 ,z N )a(X 1 ,z N ) 

N-l 

b(X N ,z N ) Y[ a(A 7 ,2:Ar) 

7=1 

(2.25) 
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It is easy to check that in the limit Ai — > A, A2 — > A, . . . , \n — » A, we have 



z N = U 



i<j 



Zj - Zj + T] 

Z{ Zj 



N 

n 



L A - Zi + 77 



a{X,Zi) ■■■ a(\,z k ) 
a N -\\,z 1 ) ••• a»-\\,z k ) 



a(X,z N ) 
a N -\X,z N ) 



(2.26) 



where and below a k (X, Zj) stands for the fc-th power of a(A, Zj). For the homogenous limit 
of the parameters Zk, we first compute the limit z 2 — > Zi = 2. Taylor expanding the second 
column as z 2 — > 2, 

G fc (A, z 2 ) = a fc (A, z) + (a fc (A, z))'(z 2 - 2) + - ^) 2 ), 

where and below, X', X" and J( n ' stand for the 1st, 2nd and 77.-th order derivatives of X 
with respect to the parameter z, respectively, and subtracting the first column from the 
second, ( I2.26P becomes 

Zi- Zj +1] ( T] \ TJ V 



( 

z N = (-v) n 



Z - Zj + 1] 



n 



Z — Zj I -*■-"■ Zi — Zi V A — Z + 77 / A A A — Z, ; + r? 



2 AT 

n 



X 



1 1 

a(A, z) (a(A,z))' a(A, Z3) 



1 

a(X,z N ) 



,r l (X,z) (a N ~ 1 (X,z))' a JV - 1 (A,z 3 ) ••• a^A^) 
Then Taylor expanding the third column as Z3 — > z, 



(2.27) 



fc (A, z 3 ) = a fc (A, z) + (a fc (A, z))' (z 3 - z) + ± (a fc (A, z))" (z 3 - z) 2 + 0((z 3 - z) 3 ), 



and subtracting multiples of previous columns, we obtain 

Z — Zj + 77 \ -p-|- Zj — Zj + 77 



n 



n 



3 AT 

n 



X 



2 j = 4 V Z ~ Z i ) jiiii Z i~ Z 3 \X-Z + T]J f^X-Zi+T] 

1 1 ••• 1 

a(X, z) (a(A,z))' (a(A,z))" a(A, Z4) ••• a(A,Z7v) 

i^-^z) (a"" 1 ^,*))' (a"" 1 ^,*))" a N -\X, Zi ) ■■■ a N -\X,z N ) 



(2.28) 



Continuing with such process, we obtain, instead of the determinant of the N x N matrix, 
the following determinant of the (N — 1) x (JV — 1) matrix for the partition function: 



n 



jV-1 
i=2 



(a(\,z))> (a(\,z))» ■■■ (a(\,z))( N -V 
(a N -\X,z))' {a N -\X,z)) n ■■■ (a N ~\X, z*))^ 



(2.29) 



To simplify the determinant in (12.291) further, we investigate its elements. Computing 



a J (A, z))^ n \ we obtain 



n k—1 



(2.30) 

k=i i=o 

where fk are functions of {(a(A, z))®} (I = 1, . . . , n) and are independent of j in (12.301) . One 
can easily obtain fk for the first and last terms, 

/ 1 = (o(A,z))W, / B = ((o(A, «))')"■ (2-31) 

Thus, by means of the properties of determinants, the partition function of the gl(l\l) vertex 
model is simplified to the following function: 

„ (-^Efel 1 * (b(\,z)) N _ 

i=2 



]JU l(a(\z)) 



(6(A,z)) 



(2.32) 



3 p/(2|l) vertex model with DW boundary condition 
3.1 Description of the model 

Let R G End(V £g> V) be the i?-matrix associated with the 3-dimensional irreducible gl{2\l) 



-Ri2(Ai, A2) 



•ading 


for V, 


i.e. 


[1] = [2] 




U3] 


= 0, 
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(3.1) 



10 



which satisfies the GYBE ( 12. ip . Here ai2, 612 and cyi are the same as those given in the 
previous section. The basis vectors |0), |1) and |2) of the 3-d gl(2\l) representation are given 
by 




Now we consider the gl(2\l) 15-vertex model on a N x N square lattice. Excluding the 
double occupancy, there are three possible electronic states, i.e. up spin f, down spin j 
and vacuum 0, at each site of the lattice. For this model, the configuration of the vertex is 
decided by the electronic states around it. Corresponding to the gl(2\l) invariance, there are 
altogether 15 possible configurations, which preserve the fermion numbers and spins, with 
non-zero Boltzmann weights. They are shown by using figures as follows. 
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Figure 3. Vertex configurations of the gl(2\l) vertex model and their Boltzmann weights. 
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Figure 4. DW boundary condition for the gl(2\l) vertex model. 



The Boltzmann weights associated with the vertices are given by the elements of the R- 
matrix (13.11) 

W\ = 1, W2 = W3 = c(A, z), W4 = . . . = Wg = a(A, z), 

w w = . . . = u>i3 = b(X,z), w 1A = w 15 = -b(X, z). (3.3) 

The gl(2\l) supersymmetric vertex model with DW boundary condition is defined as 
follows. At the left and top ends, all electrons are in vacuum states while at the right and 
bottom ends, there are P down spin states (corresponding to the i s t_pth lines) and N — P 
up spin states (corresponding to the (P + l) -iV^ n lines). The boundary condition is shown 
in figure 4. The DW partition function of the gl(2\l) vertex model is given by 

15 

i=i 

P-l P-l N-l N N P N N 

= n t(jv-i) n kp-s) n n n ^ u yi n ^« 

i=0 j=0 i=0 k=l 1=1 j=l j=P+l i=l 

P-l P-l N-l N P N N 

= n k*-i) n kp-j) n v^n^^n ^ n 

i=0 j=0 i=0 k=l j=l j=P+l i=l 

12 



P-l P-l N 

= n t(*-o n i(p-,)^(Ai)...c 2 (A P )c 1 (A P+1 )...c 1 (A iV )n^) ) 

i=0 J=0 8=1 

(3.4) 

where are the number of configurations with the weights Wi, and T k (X k ) is the monodromy 
matrix along the horizontal lines and is defined by 

T k (X k ) = i?fc,jv(Afc, 2;jv)-Rfe,jv-i(Afc, zn-i) ■ ■ ■ Rk,i(Xk, z\) 
An(X)k) A 12 (X k ) ^(Afc) \ 

A 21 (A fe ) A 22 (A fe ) F 2 (A fe ) . (3.5) 
CiCA*) C 2 (A fc ) D(A fc ) / (fe) 

Similar to the pZ(l 1 1) vertex model case, the partition function (13.41) can be computed by 
using the approach of Drinfeld twists, as can be seen from the next subsection. 

3.2 Exact solution of the partition function 

We now compute the partition function H3.4[) using the Drinfeld twist method. The F-matrix 
for this case is still defined by ( 12.101) . except that now a = 1, 2, 3. The inverse of the F-matrix 
is given by 

i<j 

with 

^ij diag (A&ijdji^ Cljij CLjij Giij, A(Xij(Xji^ Clji^ (^ij^ CLij , 1) - (3.T) 

Working in the F-basis, the lower entries C\ and C 2 of the monodromy matrix (13.51) are 
simplified to symmetry forms, that is they can be written as [19] 



N 

C 2 (X) = ^&(A,^)F ( 2 ; } 3 diag(a(A,Zj),2a(A,Zj)> l) (j) , 
i=i 

(3.8) 

N 

= ^2 b (Kzi)E((l diag (2a(X, zj), a(X, Zj)(a(zi, Zj))~ l , l) (j) 

8=1 



V> a(A, Zj)b(X, Zj)b(zj, zj) 1>2 f2;3 
^ a(z h z 3 ) E ® E ® 

®k^i,j (2a(X,z k ),a(X,z k )a' 1 (z i ,z k ), l) (fc) . (3.9) 
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Here (a, {3 = 1, 2, 3; I — 1, 2, • • • , N) are the generators of gl(2\l) on the Zth vertical line. 

In the following we will identify |0), |1) and |2) with the vacuum <f), spin-up f and spin- 
down I states, respectively. Note that under the action of F the state |0)(i) ® • • • ® \0)(n) is 
invariant, that is, 

Fi...7v|0) ( i) <g> • • • ® |0) (JV ) = |0) ( i) ® • • • <g> |0) (JV ). (3.10) 
Thus substituting fl3.8p - fl3.10p into the partition function (13.41) . we obtain 

P-l P-l JV 

= n( i i(^)n( 2 i(^-^(Ai)...c 2 (Ap)c 1 (A P+1 )...c 1 (A^)nio) W 

i=0 j'=0 i=l 

P-l P-l 

= n^i(^) n^i^^-^cAi) . ..c 2 {\p)c 1 {x P+1 ) . ..cm n io> (0 

1=0 j=0 1=1 

P-l P-l P N N 

= e e n<v-on<%-i) F ^n^ n^iii% 

ii<...<ip ip + i<...<iAT 2=0 i=0 j=l j=P+l 1=0 

P N 

P(P+1) (JV-P)(JV-P+1) -p-p -p-,- 

fc=l Z=P+1 

xdetB P ({Ai, • • • , Ap}, {z h . . . , z ip }) 
xdetBjv_p({Ap + i, . . . , Ajv}, {zi P+1 . . . , z iN }) 

P N 

= E E XiiM^^r 1 II (2a(^,^))- 1 G({ 2 }) 

il<...<ip i P+ i<...<ijv l>k=l l>k=P+l 

P N 

P(P+1) { N-P)(N-P+1) yj yj 

x2 2 + 2 II II "(A/,- "-J 

fc=i z=p+i 

xdetBp({Ai, • • • , A P }, . . . , z ip }) 

xdetBjv-p({Ap + i, . . . , Atv}, \Zi p+1 . . . , z iN }), (3.11) 

where {ii, . . . , ip} n {ip+i, . . . , i n } = and £>m({A}, {z}) is a M x M matrix with elements 

«-i 

0M{A}, {2}))a^ = 6(A a , ^) n «(A 7) */»)• (3-12) 

7=1 

In ( 13. lip , the function G({z}) is defined by 

P N 

G({z}) = l[(2a(z lk ,z H )) J] (2a(^,^))(l| (A r)...(l| ( p +1) (2| ( p)...(2| (1) 

£>A;=1 2>fc=P+l 

• • • E? p El? p+i) ,E\! N) |0) (1) . . . |0) (AO (3.13) 
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Using the matrix F 1 ] N defined by (I2.15p - fl2.18p . one may simplify G({z\) as follows. 
p N 



G({z}) = Y[(2a(z tk , Zil )) J] (2a(z ik , z k )) 

>k=l l>k=P+l 

: (0| (JV ) • • • (0| ( i)^) • • • E th + i) E (P) . . . Eg S(c, a, a a ) 



l>k=l [>k=P+l 
X \ 

^iVwl?) • • • Km II A ^ 4o • • • • • • 10) w 

i<j 

P N P N 

= n n a_i (^>^) n cm^.*.)) -1 n (M^^r 1 

fc=lZ=P+l fc>Z=l fc>Z=P+l 

x (0| (JV) . . . (OldjE^, . . . ^ +1) E ( 3 |) • • • Eg E 5 ^ c ' a ' 

<tS<Sjv 

X K(L.N) E fh) ■ ■ ■ E ip E (i P+1 ).. E (iN) ■ ■ • 

P N P N 

= n n a-i (^fc>^) n (m^^)) -1 n ck**^)) -1 

A;=1Z=P+1 fc>/=l fc>Z=P+l 

x e (-d^w «-o «»s(c,., Q ,) {ki^)2~°" • (3l4) 

where a = 1, 2, the subscribes of a are indices of space, sign(a) = 1 if a is odd and sign(a) = 
if a is even. 

The procedure of computing the homogenous limit is similar to that for the gl(l\l) vertex 
model. Here we only give the results. In the homogenous limit, i.e. when Ai = . . . = A^v = A 
and Z\ = . . . — zn = z, the partition function (13. 1 1H becomes 

Z N = (a(X,z)r N - p \b(X,z)r(b(X,z)Y N -^ 

x lim E G ^ ( 3 - 15 ) 

ii<...<ip ip+i<...<ijv 

By using the decomposition law (12.131) . the .R-matrix R a in (I3.14p can be decomposed to 
elementary R- matrices. In homogenous limit, the elements of the elementary i?-matrix will 
tend to 

R(z, z)H = R(z } z)H = -1, R(z, z)\\ = R(z } z)\% = -1, R(z, z)\ 2 2 = R(z } z)\\ = 0. (3.16) 
Therefore for the last factor in (13.151) involving G({z}), one may easily obtain 



lim E E G ^ 

«i<...<ip «p+i<...<ijv 
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P N 




z yi e ini fl \ z i k > z h) 



il<...<ip ip+i<...<ijv k=l l=P+l 




(3.17) 



where is the combinatorial number. 

In summary, in the homogeneous limit the DW partition function Zn is 



(3.18) 



It is easily seen that when P = or P = N, the DW partition function of the gl(2\l) 
vertex model reduces to that of the g/(l|l) vertex model, confirming a statement made in 
[21] on DW partition functions of vertex models based on A^-type algebras. However, when 
P 7^ 0,N, one obtains new partition functions of the true g/(2|l) supersymmetric vertex 
model with DW boundary conditions. 

4 Conclusion and discussion 

In this paper, we have proposed the gl(l\l) and gl(2\l) supersymmetric vertex models with 
the so-called DW boundary conditions. The DW partition functions of the models have 
been computed by means of the approach of the Drinfeld twists. We have found that in 
the homogenous limit, the partition functions degenerate to simple functions. For the 2-d 
square gl(2\l) lattice model, we note here that the definition of the DW boundary condition 
is not unique. For the other cases, by using the same procedure, one may find their partition 
functions are similar with those in this paper. 

We have demonstrated, by working out the DW gl(l\l) and gl(2\l) supersymmetric vertex 
models as examples, that by means of the Drinfeld twist method, one can actually derive 
directly, rather than conjecture a formula and then verify it as usually done, the determinant 
representations of DW partition functions. 

It is widely known that determinant representations of partition functions are closely 
related to some pure mathematical problems, such as algebraic combinations and tilings of 
the Aztec diamond [51 [22]. In our further work, we will study the mathematical problems 
arising from the present models. The results in this paper will also be useful for simplifying 
the correlation functions of the supersymmetric t-J model obtained in [23J and for further 
studying physical properties of the model. 
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